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The possibility of novel behavior at interfaces between strongly and weakly correlated materials 
has come under increased study recently. In this paper, we use determinant Quantum Monte Carlo to 
determine the inter-penetration of metallic and Mott insulator physics across an interface in the two 
dimensional Hubbard Hamiltonian. We quantify the behavior of the density of states at the Fermi 
level and the short and long range antiferromagnetism as functions of the distance from the interface 
and with different interaction strength, temperature and hopping across the interface. Induced 
metallic behavior into the insulator is evident over several lattice spacings, whereas antiferromagnetic 
correlations remain small on the metallic side. At large interface hopping, singlets form between the 
two boundary layers, shielding the two systems from each other. 

PACS numbers: 71.10.Fd, 71.30.+h, 02.70.Uu 
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I. INTRODUCTION 

Over the last decade, advances in synthesis techniques 
have made possible the construction of well-defined in- 
terfaces involving strongly correlated materials, notably 
transition metal oxides^—. Indeed, the early discovery 
that the interface between two insulating oxides, LaA10*3 
and SrTiC>3 is a high-mobility two dimensional conduc- 
tor— or even a superconductor— has emphasized that novel 
physics can arise at such boundaries, beyond a simple in- 
terpolation between materials properties on cither side. 
Subsequently, numerous heterostructurc interfaces have 
been shown to exhibit unique phenomena that are not 
present in their bulk constituents. In particular, Mu- 
nakata et al. presented a subtle proximity effect that 
arises between a normal metal and an antiferromagnetic 
Mott insulator, which can be understood in the frame- 
work of Rudcrman-Kittcl-Kasuya-Yoshida (RKKY) in- 
teractions 3 -. 

Advances in characterization techniques are also ex- 
panding the horizon of the field. It was shown recently 
that, by setting up a standing wave of the incident X- 
rays and adjusting the position where the intensity is 
peaked inside the sample, one can measure electronic ex- 
citations as a function of location in the sample. This 
standing wave angle resolved photoemission spectroscopy 
(SW-ARPES)£ method allows a depth selective probe of 
buried layers and interfaces, and hence the construction 
of chemical and electronic structure profiles within the 
sample. 

Various numerical methods have been employed to 
model qualitative features of these experiments. Early 
attempts include the study of LaA10*3 and SrTiO*3 
systems in which the electric fields arising from the 
La 3+ /Sr 2+ charge difference and the on-site interactions 
are treated within the Hartree-Fock approximation^. A 
metallic interface between band and Mott insulators in 
a quasi-one-dimensional lattice was explored with the 
Lanczos method^. Important insights into electronic 



properties at interfaces have also been gained by dynam- 
ical mean field theory (DMFT) 1 ^. For example, it has 
been suggested that the Kondo effect governing the in- 
terface between metal and Mott insulator is inefficient 
so that Mott insulators are impenetrable to the metalii. 
More specifically, the quasiparticle weight decays as 1 jx 1 
with distance x from the metal, but the prefactor of this 
decay was found to be very small. Zenia et al., however, 
emphasized that even a small proximity effect can induce 
density of states to open up a metallic channel inside an 
insulator sandwiched between to metals at sufficiently 
low temperatures 1 ^, leading to perfect conductance. This 
Fermi liquid is "fragile" : finite temperature, disorder, or 
frequency rapidly return the behavior to that of a conven- 
tional N-I-N junction. Possible device applications were 
suggested as a consequence of this sensitivity^— In a re- 
lated worki^, a Gutzwiller approximation approach was 
extended to inhomogeneous systems. The decay length 
of the exponential fall-off of the penetration of metallic 
character into the insulator region was shown to diverge 
as the metal-insulator transition is approached. 

A further interesting set of studies involves the dy- 
namical response of strongly correlated systems with 
an interface, for example the possibility of the sup- 
pression of charge transport, current rectification, and 
the behavior of holon-doublon pairs, which are impor- 
tant for devices^. So far, such dynamic phenomena 
have been explored primarily for one-dimensional sys- 
tems, where time-dependent Hartree Focbi^ and density- 
matrix renormalization group (DMRG) methods are es- 
pecially effective^. 

In this paper, we address the metal-insulator interface 
problem by using determinant Quantum Monte Carlo 
(DQMC) 17 ' 18 , a numerically exact approach, to treat 
correlated electron models. Here we focus on an inho- 
mogeneous Hubbard model on a single two dimensional 
lattice divided by a linear interface in two regions, one 
weakly correlated (with U = 0), and one at intermedi- 
ate coupling (with U non-zero). The chemical potential 
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is chosen to maintain particle-hole symmetry, so that all 
sites of the lattice are on average half-filled. This choice 
avoids the fermion sign problem^ and enables the eval- 
uation of magnetic correlations and the density of states 
at low temperature. 

Such DQMC simulations are at present limited to lat- 
tices of 400 — 1000 sites (depending on the interaction 
strength and temperature). By focussing on a linear in- 
terface in a 2D lattice, we are able to explore systems with 
a fairly large linear extent. Most of our results will be for 
20 x 20 lattices. This enables us, for example, to evalu- 
ate the penetration depth across the boundary, since the 
interface effects have sufficient room to heal before the 
lattice edge. The behavior of the Hubbard model across 
a planar interface in a 3D material on systems of smaller 
linear extent has been explored by Euverte et a&^. The 
key finding is that up to an interface hopping V which is 
on the order of the bulk hybridization t, the effects of the 
interface can extend well past the two layers immediately 
at the interface. That is, the interface affects properties 
3 — 4 layers deep on both the strongly and weakly cor- 
related sides of the boundary. When V > 2t, however, 
there is a return to the values characteristic of decoupled 
materials with V = 0. This "revival" of magnetic or- 
der on the insulating side is driven by the formation of 
singlets between fermions on the two layers immediately 
adjacent to the interface which acts to decouple the two 
materials. 

This paper is organized as follows: In section II we ex- 
plicitly write down the Hamiltonian and provide a very 
brief overview of DQMC. Section III focuses on the mag- 
netic properties of the 2D lattice with a linear interface, 
both the near-neighbor spin correlation and the structure 
factor. Here the interaction strength U and temperature 
T are varied at fixed interface hopping V. The key result 
is a determination of the penetration depth as a function 
of T and U. Section IV extends these results to the den- 
sity of states. Section V examines the variation with V, 
and Section VI summarizes our results. 



II. MODEL AND METHODOLOGY 

We consider the two-dimensional Hubbard model 

H = -t ^ ( C iL<T C jL<T + hM -) - M niL ° 
(ij) ,L,cr i,L,cr 

- t LL'(4 La C lL , lJ + h.C.) 

i,(LL'),a 



(i) 



i,L 



Here c\ Lcr (c iLcr ) are fermionic creation (destruction) 
operators of spin a at site i in line L. While i, L can 
also be regarded as the x and y site labels i x ,iy, the cur- 
rent notation emphasizes the broken rotational symmetry 
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FIG. 1: (color online) Geometry of the two-dimensional lat- 
tice with a one-dimensional interface (dashed red box). The 
metallic lines (negative L) have U — while the insulating 
lines (positive L) have U 7^ 0. The lattice has periodic bound- 
ary conditions in the x direction and open boundary condi- 
tions in the L direction. V denotes the hopping/hybridization 
(heavy lines) across the interface. We adopt the lattice size 
N = 20 x 20 throughout the paper, unless otherwise indicated. 



(the interface is parallel to the x axis) . Our labeling con- 
vention is L = ■ ■ ■ — 3, —2, — f for the U = lines and 
L = 1, 2, 3, • • ■ for the U 7^ lines, so that the interface 
connects L = ±1. t and tw are the intra and inter-line 
nearest-neighbor hoppings. tw = t except at the inter- 
face where t-\\ = V . The geometry is shown in Fig. 

ED 

The DQMC approach provides an exact solution to 
the finite temperature properties of H by decoupling the 
interaction term U through the introduction of an aux- 
iliary Hubbard-Stratonovich (HS) field. The resulting 
action is quadratic in the fermion operators, so that the 
trace over those coordinates can be performed, leaving 
determinants (one for spin up and one for spin down) 
of matrices whose dimensions are the spatial lattice size, 
and whose entries depend on the HS field. The auxil- 
iary field is sampled stochastically, and the values of the 
up and down spin fermion Green's function for the con- 
figurations generated are used to construct the various 
observables. The CPU time for the algorithm scales as 
the cube of the lattice size and linearly^ with inverse 
temperature /3. 

At low temperatures, the fermion determinants can be- 
come negative, precluding their use as a probability for 
sampling the HS field. To avoid this "sign problem" , 
we consider the half-filled case, p, = 0. Here particle- 
hole symmetry (PHS) implies that (riiLo) = 0.5 for all 
temperatures and interaction strengths, even if Ul vary 
spatially with L and the hybridizations are not all equal. 
PHS can also be used to demonstrate that the up and 
down spin determinants, although they can individu- 
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FIG. 2: (color online) Spin correlation Cl{x) = {Si,lSi+x,l) 
as a function of separation x for different lines L. Diminishing 
antiferromagnetism of the insulator adjacent to the interface 
is shown. The induced antiferromagnetic long-range orders in 
the metallic region are evident in spite of their weaknesses. 
U = 0|4 denotes the interface setup. 



ally become negative, always have the same sign. As 
a consequence, their product is always positive, allowing 
the study of the low temperature behavior at the metal- 
insulator interface. 



III. MAGNETIC PROPERTIES 

We gain quantitative information concerning the spin 
correlations in different lines by measuring the spin corre- 
lation C L {x) = (S~ L S+ L ), with Sr L = c n 4 p parallel 
to the interface. Figure [3] shows that antiferromagnetism 
in the insulating layer L = 1 immediately adjacent to 
the metal is diminished by 20-30% relative to larger pos- 
itive L which are deeper in the Mott insulator. Data are 
shown for boundary hybridization V — t and interaction 
U = At in the insulator. Despite a somewhat smaller am- 
plitude, the correlations in L = 1 appear to remain long 
ranged. Induced antiferromagnetic long-range order— in 
the metallic region is also evident, although its ampli- 
tude is an order of magnitude smaller than that on the 
insulating side. 

Compact information on the antiferromagnetic long- 
range order for each line L can be obtained by measuring 
the antiferromagnetic (AF) structure factor, 

SAF{L) = ± y £(-ir-C L (x) (2) 

x 

where N is the linear lattice size. Figure [3] illustrates the 
L-dependent AF structure factor for different correlation 
strengths U and temperatures T. The smooth penetra- 
tion of antiferromagnetic order into the metal and the 
diminishing AF order in the insulator are clearly seen. 
As might be expected, the induced long range AF order 
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FIG. 3: (color online) AF structure factor Saf(L) in each 
line L for different correlation strengths U and temperatures 
T. The smooth penetration of antiferromagnetic order into 
the metal and the diminishing AF order in the insulator 
are clearly seen. Fitting Saf{L) to a hyperbolic function^, 
atanhi/A + b allows the extraction of a penetration length 
A, shown in the insets. A decreases with larger correlation U 
and lower temperature T. Further information 



in the metal adjacent to the interface is stronger for corre- 
lated insulator with larger— on-site repulsion U (Fig. [3^) 
and/or lower temperature (Fig. [3}d). The overall magni- 
tude of Saf in the metal is rather small, consistent with 
the small real-space correlations in Fig. [3J The results of 
Figs. 12131 arc rather different from previous work of Sher- 
man et. al^ which found that the antiferromagnetic 
order can penetrate into the metal to a depth of ten lat- 
tice spacings, but are consistent with the shorter range 
effects described i n 11 ' 12 . Figure [3] indicates that by layer 
L = — 3 the influence of the contact with the insulator is 
minimal. 

On the other hand, we find that contact with the metal 
has a substantially greater effect on the insulator. Here 
the "Kondo proximity effect" diminishes the long range 
AF order of the insulating lines adjacent to the inter- 
face. This tendency to paramagnetism competes with 
the long range order induced by the exchange energy 
J ~ t 2 /U deep in the insulator (Fig. [3]i). To quantify 
the proximity effects, wc fit the curves of Saf{L,U,T) 
with the hyperbolic form Saf(L, U,T) = a tanh L/ X + b. 
The penetration depth A decreases with larger corrcla- 
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FIG. 4: (color online) The L dependence of the density of 
states N(uS) shows the evolution of the spectral properties 
across the metal-insulator interface. Lines of sites within the 
Mott insulator (L > have nonzero U) are characterized by 
the presence of Kondo resonance peaks split by antiferromag- 
netic order at half-filling. However this AF gap between the 
two peaks is dramatically weakened as L —} 1. On the other 
hand, the densities of states associated with lines of sites on 
the metallic side of the interface show little effect of contact 
with the insulator. Even for the line L = — 1 most immedi- 
ately adjacent to the boundary there is only a small dimple 
in N(oj) at uj — 0. 

tion U. Interestingly, A increases at higher temperature. 
It is possible that raising T weakens the magnetic or- 
der in the insulator, thereby enhancing the effects of the 
contact with the metal. 



IV. DENSITY OF STATES 

One criterion to distinguish metal from insulator is the 
single particle density of states, which can be extracted 
from the local imaginary-time dependent Green's func- 
tion G l (t) =E^x.W4 ff (0) from 

to* 

/OO p-CJT 
_ d --^TT Nl(u3) (3) 

and using the maximum entropy method^. The density 
of states is probed experimentally with photocmission 
spectroscopy. In a translationally invariant system one 
often also calculates the momentum dependent spectral 
function A(k,u), measured in angle-resolved photoemis- 
sion spectroscopy (ARPES). 

As mentioned in the introduction, a dramatic experi- 
mental achievement in the past few years has been the 
development of standing wave ARPESi which has en- 
abled the probing of strong correlation effects layer by 
layer in a sample. Here we make some predictions for 
possible features to be seen in SWARPES by obtaining 
the row dependent single particle density of states. As 
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FIG. 5: (color online) L-dependent density of states for dif- 
ferent on-site repulsions U and temperatures T. The smooth 
penetration of metalicity into the insulator and diminished 
metallic behaviors in the metal are clearly seen. Similar to 
Fig. O we employ a hyperbolic fitting function and extract 
a penetration depth A (insets), which we find decreases with 
larger correlation U and increases with higher temperature T. 



with the row- by- row magnetic correlations, Nl(lo) shown 
in Fig. [4] characterizes the penetration of metalicity into 
the insulator. Kondo proximity effects are evident in the 
evolution of Nl(oj). Figure |4] shows that the insulating 
lines for L = 6, 5,4, relatively far from the interface, are 
characterized by the presence of a Kondo resonance peak 
split by antifcrromagnctic order at half-filling. However 
as the interface is approached for L = 3, 2, 1 the gap 
is increasingly filled in, evidence of the coupling to the 
metallic half of the lattice. Meanwhile the metallic line 
[L = —1) immediately adjacent to the interface shows 
some influence of the boundary, albeit a rather small one: 
the central peak at uj — shows a slight dip. This is com- 
pletely gone for L = — 2. Evidcntally the metallic behav- 
ior of N(ui) penetrates much further into the insulator 
than the insulating physics does into the metal. 

Further information on the spectral weight at Fermi 
level Nl(U, T)\ u —o as a function of on-site repulsion U 
and temperature T is shown in Fig. [5] The penetration 
of metallic behavior into the insulator and, conversely, 
the diminished spectral weight at the Fermi surface in 
the metal are clearly seen. Similar to Fig. [31 wc adopt a 
hyperbolic function for fitting. The resulting penetration 
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FIG. 6: (color online) (a) Evolution of the antiferromagnetic 
(AF) structure factors Saf{L) with interface hybridization 
V. Except the line of sites L = 1, adjacent to the metal, 
Saf(L > 1) exhibits a nontrivial revival with increasing the 
V. We attribute this affect to singlet formation at the metal- 
insulator interface which leaves the lines L > 1 decoupled 
from the metal, (b) Further evidence of the strong magnetic 
coupling across the interface. U = 0|4, /3 — 10. 



depths A are given in the inset, and decrease with larger 
correlation U and lower temperature T. 



V. EFFECT OF VARIATION OF INTERFACE 
HOPPING 

We show in Fig. [5] the evolution of the line dependent 
AF structure factors Saf(L) with interface hopping V. 
As expected, all lines L > with nonzero on-site repul- 
sion U = 4 have large Saf when the metal and insulator 
are decoupled (V" = 0). With increasing hybridization all 
Saf{L > 0) initially decrease. The most dramatic fea- 
ture in Fig. inja) is the different behaviors of Saf{L = 1) 
and Saf(L > 1). In the line L = 1, adjacent to the 
metal-insulator interface, Saf{L = 1) monotonically de- 
creases with V. However, Saf{L > 1) show nontrivial 
revivals with increasing V and ultimately recover to val- 
ues characteristic of the decoupled case V = when there 
is no contact with the metal whatsoever. This behav- 
ior of Saf{L) is qualitatively similar to the multi- layer 
metal-insulator interface in the three-dimensional Hub- 



bard models, and can be attributed to "singlet" for- 
mation at the metal-insulator interface, which both sup- 
presses Saf{L = 1) and also leaves the remaining lines 
L > 1 decoupled from the metal. 

The "metallic" line structure factors Saf(L) retain 
their small values for the entire range of V. As with the 
real space spin correlation functions, there is a greater 
effect of the metal on the insulator (suppressing mag- 
netic order) than of the insulator on the metal (inducing 
magnetic order). 

Figure [S^b) provides further evidence of the strong 
magnetic coupling across the interface. The dominant 
feature is the rapid increase of antiferromagnetic cor- 
relation across the metal- insulator interface C±(L = 
— 1, V = 1) as increasing the hopping V. Note also that 
magnetic correlation C±(L = 1,L' = 2) become smaller 
with larger hopping V due to the "singlet" formation at 
the interface L = —1,1. 



VI. CONCLUSION 

We have used the numerically exact finite-temperature 
determinant Quantum Monte Carlo (DQMC) method to 
study a metal-insulator interface in a two-dimensional 
square lattice. By investigating the long-range antiferro- 
magnetic order and density of states, we demonstrated 
that the metallic behavior penetrates into the insulator 
for several lattice spacings, with a penetration depth A 
which decreases with increasing on-site repulsion, U, on 
the insulator side of the interface. A is also (somewhat 
more weakly) temperature dependent, decreasing as T 
is lowered towards the critical temperature T c = for 
long range magnetic order on the insulating side. The 
penetration length A thus shows an opposite tempera- 
ture dependence to that of the spin correlation length, £, 
which grows as T is lowered. 

The insulator-induced antiferromagnetic long-range 
order in the metal is predominantly limited to the line 
immediately adjacent to the interface. That is, magnetic 
characteristics of strong correlation appear to penetrate 
less deeply into the metallic side of an interface than 
does weak correlation physics penetrate into the insula- 
tor. This is consistent with previous results for a planar 
interface in a 3D lattice^. We note, however, that in the 
3D geometry, the effect of the contact with the insula- 
tor on the in-plane conductivity a in the U = half of 
the lattice can extend beyond the contact layer. This is 
consistent with our results for the density of states at the 
Fermi level N(ui = 0), which arc shifted from their U = 
values for layer L = — 2 in addition to layer L = — 1. 

In the past several years, studies of the effect of 
spatially varying densities and interaction strengths 
have been motivated by experiments on trapped atomic 
gase&2£~—. In such systems the spatial variation, eg a 
quadratic confining potential Vt rap (r) = Vr(r/l) 2 , has 
an explicit length scale I. This complicates the deter- 
mination of intrinsic length scales associated with the 
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response of the interacting fermions themselves. Our 
choice here of a sharp (scale free) interface between 
metal and insulator (Ul — 0, L = —1,-2,-3... and 
Ul = U > 0, L = 1,2,3...) allows us to attribute the 
lengths characterizing the relaxation of properties on ei- 
ther side of the interface solely to the fermionic correla- 
tions. This choice is, of course, also more appropriate to 
the solid state context of a sharp interface between two 
materials. 

One further motivation to study the metal-insulator 
interface is its close relation to the question of "orbitally 
selective Mott transitions" (OSMT)2Ir^&. Here the cen- 
tral question is whether orbitals with different degree of 
electronic correlation, coupled together by interorbital 
hybridization, necessarily undergo the Mott transition 
simultaneously. The layer index 'L' in the Hamilto- 
nian considered in this paper bears a formal similarity 
to the orbital degree of freedom in the OSMT, although 
of course the details of the coupling via hybridization 
are rather different in the two cases. It is evident from 
our data that a layer which shows the hallmarks of an 
antifcrromagnetic Mott insulator can coexist with layers 
which have the characteristic behavior of a paramagnetic 
metal. That such coexistence is possible is similar to the 



conclusion ultimately reached in numerical studies of the 
OSMT. 

In this paper, the density of states and magnetic cor- 
relations near an interface between {7 = and U ^ 
regions have been found to be more or less smooth inter- 
polations between the bulk metal and Mott insulator. On 
the other hand, Quantum Monte Carlo simulations have 
observed novel phases^ such as spin liquids, to arise in 
models where the energy scales are poised at the bound- 
ary between a semi-metal and an antiferromagnet. An 
interesting extension of the present work would be, there- 
fore, to bring such more general regions into contact and 
study the properties at the interface. 
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